Deceptive signals of phase transitions in small magnetic clusters 
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We present an analysis of the thermodynamic properties of small transition metal clusters and show how 
the commonly used indicators of phase transitions like peaks in the specific heat or magnetic susceptibility 
can lead to deceptive interpretations of the underlying physics. The analysis of the distribution of zeros of 
the canonical partition function in the whole complex temperature plane reveals the nature of the transition. 
We show that signals in the magnetic susceptibility at positive temperatures have their origin at zeros lying at 
negative temperatures. 

PACS numbers: 64.60.-i, 36.40.Ei, 05.70.Fh 



o 

e 

I 

-(— > 



X3 

o 

> 

o 



cd 



I 

O 

o 



X 



Experiments on Bose-Einstein condensation [111 0] or the 
experimental determination of structural, electronic and ther- 
mal properties of clusters [Q §, ||| are prototypes of physi- 
cal investigations of transitions in small systems. Intuitively 
phase transitions in such systems do exist. While the atomic 
structure at low temperatures is more or less rigid, at high 
temperatures the atoms move resembling a liquid drop. But 
the theoretical description is complicated since the thermody- 
namic functions of clusters do not show singularities at the 
transition point. Phase changes are seen in blurred slopes or 
humps. To have the physical concept of phase transitions and 
their properties of bulk material in mind and apply it to inter- 
pret the smooth thermodynamic functions for small systems, 
e.g. clusters, of a given size accordingly and assign a "first" or 
"second" order phase transition may be inconclusive. Despite 
the ambiguity in these signals, it is still reasonable to attribute 
an order to phase changes in small systems because funda- 
mental differences between the kind of the transitions such as 
the existence of metastabilities for first-order transitions still 
persist. Therefore, various approaches for the classification of 
phase transitions in small systems have been developed which 
have to coincide in the thermodynamic limit and should be 
mathematically rigourous. 

Gross et al. have proposed a microcanonical treatment, 
where phase transitions are distinguished by the curvature of 
the entropy S{E) [0, ^. If S{E) has a convex intruder, i.e. the 
microcanonical caloric curve T{E) shows a backbending, the 
phase transition is assumed to be of first order Franzosi et 
al. have started by investigating the topology of the potential 
energy surface and established a connection between topolog- 
ical changes and phase transitions [& |lO|. However, they are 
not able to determine the order of the phase transition. Re- 
cently we have proposed a classification scheme based on the 
distribution of zeros of the analytically continued canonical 
partition function Z{B), with B = f3 + ir (j3 = l/T), in the 
complex temperature plane [ [ll| , [I4 , [l3[ [l4| ] . 

The basic principle of the description of phase transitions 
by the zeros of the partition function is the product theorem of 
Weierstrass and the theorem of Mittag-Lefiffer which relate in- 
tegral functions to their zeros Jl^. Applying these theorems. 



the canonical partition function can be written as 
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We assume the zeros to lie on a line with a density (/)(r) ~ 
r" and to have a crossing angle v with the plummet on the real 
temperature axis (7 — tan v). Together with the imaginary 
part Ti of the first zero Bi this leads to a distinct classification 
of phase transitions in small systems. For zeros perpendicular 
to the real axis with equal or increasing spacing, i.e., a < 
and 7 = 0, the transition is of first order, for < a < 1 and 
arbitrary 7 as well as for 7 7^ and a < of second order, see 
Fig. nl The imaginary part ti reflects the "discreteness" of the 
system. Thus, in the thermodynamic limit we have ti ^ 
and our scheme coincides with the scheme given by Gross- 
mann and coworkers |p^] . 
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FIG. 1: Examples of distributions of zeros for P' and 2"^* order 
phase transitions along with the corresponding specific heat, calcu- 
lated as functions of the zeros. 



We utilize small magnetic clusters exposed to an external 
magnetic field in order to show how the common treatment 
of phase transitions in small systems like the identification by 
humps of response functions eventually leads to misinterpre- 
tations of physical properties. 

Metal clusters have the intriguing property that they occur 
as different isomers with almost equal ground state energies 



but very different magnetic moments and different geome- 
tries [|l7| |l8||. For simplicity, we consider in our model two 
isomers with magnetic moments /ii = I/ib and /i2 = lOyUe. 
and their ground state energy difference AE = Eq{2) — 
-Eo(l) = 1 meV. In the presence of an external magnetic field 
H pointing in z-direction the partition function reads 
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We have assumed equal vibrational energies. The two isomers 
can be identified by their average magnetic moment (/i) which 
are calculated by standard differentiation of Eq.(||) with regard 
to the magnetic field 

(/i) = /3-1 dn lnZ(/3) = J2p^ {^i^) (4) 

i 

- Y,P^ bHtanh-\f3fi,H) - l/{(3H)l (5) 

i 

where pi is the probabihty of finding isomer i. 
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FIG. 2: Contour plot af the average magnetic moment (/i) versus 
temperature T and magnetic field H . Note the non-linear scale. 

This system is driven by two effects, the entropy increase 
due to thermal excitation and the alignment of the magnetic 
moments along the magnetic field. For low temperatures, 
there is a transition from (/i) ~ 1 to (/i) ~ 10 at about 
1.1 mGauss, as shown in Fig. g. At higher temperatures the 
magnetic field does not align the magnetic moments along the 
field resulting in a general decrease of (/i) . 

Figure |3Ja) shows the occupation probability of isomer 1 
(/ii = I^b) for different magnetic fields. At low magnetic 
fields, the lower energetic isomer is predominantly occupied, 
while higher magnetic fields lower the ground state energy of 
isomer 2 (/i2 = lO/is) and therefore the occupation is re- 
versed. With increasing temperature the probabilities become 
equal. The contributions of both isomers to the total average 
magnetic moment (/i) are plotted in FigB(b) and (c). At low 
temperatures, small magnetic fields align the magnetic mo- 
ment of isomer 1 H. With increasing temperature the mobil- 
ity of the atoms is raised which decreases the contribution of 
isomer 1 to (/i). 



From Fig. || and Fig. |] we can infer that for temperatures 
T < 1 K and an increasing magnetic field a transition with 
a coexistence phase occurs (for coexistence in small systems 
see [p^pC, 21 , 22[|). As the magnetic field is raised the con- 
tributions of both isomers to (/i) are inverted which causes a 
bimodal probability distribution P{fiz) of the order parameter 
fiz- The response of the system is observable as a hump in the 
susceptibility x = dnifJ-) at about 30 K for H = 2.0 mGauss, 
see Fig. |(d). 
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FIG. 3: Plots of physical properties versus temperature T for H — 
0.1, 0.5, 1.0, and 2.0 mGauss. (a) Probability for isomer 1. (b) 
Contribution of isomer 1 and (c) of isomer 2 to the total average 
magnetic moment (/i) . (d) Magnetic susceptibility x- 

However, at temperatures about 10 — 1000 K the situation 
is a bit more complicated. With increasing temperature and 
at "intermediate" fields (H < 1.0 mGauss) the contribution 
of isomer 1 decreases, while the contribution of isomer 2 to 
(/i) for is raised up to 1 /i^. This also results in humps of the 
susceptibility x but P{fJ.z) is unimodal because the magnetic 




p[1/eV] 



p[1/eV] 



p[1/eV] 



\^>mB] 



a 



FIG. 4; The absolute value of the average magnetic moment | {fi) \ in the complex temperature plane for (a) O.lmGauss, (b) 1.2 mOauss, and 
(c) 2.0 mGauss. 
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FIG. 5: Contour plot of the magnetic susceptibility x over the complex temperature plane, (a) for H = 0.5 mGauss, (b) displays a close-up of 
ca. one tenth of this plane, and (c) for H — 2.0 mGauss. The white crosses indicate the location of the poles of x- The arrows are a guide to 
the eye and visualize the "radiation" of the poles. 



moments of the lO/is isomer are less aligned along the mag- 
netic field. Since the first effect can be regarded as pure mag- 
netic field driven, it is not clear whether one should attribute 
the humps for H < 1.0 mGauss to a magnetic field effect or 
to a temperature effect. The "coexistence" observable in the 
contributions of both isomers to (//) is fundamentally differ- 
ent, it would be more coiTect to associate this with thermal 
excitation. 

Magnetic clusters are finite spin-systems. Such systems, if 
their energy has an upper limit, can show an inverse change 
of entropy with respect to energy coiTesponding to negative 
inverse temperatures, (3 = 1/T = dES{E) = 9b In ri(i?) 
with il{E) being the density of states. Note, that negative 
temperatures are attributed to the spin temperature of the sys- 
tem. Negative temperatures have been measured, e.g. in LiF- 
crystals [ |23| , p4| , |2^ by decoupling the spin temperature from 
the kinetic energy contribution. The spin-lattice relaxation 
time is large enough (up to hours) to assure that the spin sys- 
tem is thermally stable and thus can come to equilibrium pq]. 

Obviously, the above presented indicators of phase transi- 
tion thwart the classification and to some extent the distinction 



between different phase transitions. Within the microcanon- 
ical ensemble the occurence of negative temperatures arises 
naturally because T is an internal parameter in contrast to the 
canonical ensemble. We consider the positive and negative 
inverse temperatures within the canonical treatment to assure 
having sufficient information. 

Figure Q shows 3-dimensional images of the distribution of 
zeros of Z{B) in the complex temperature plane. The poles of 
{^) I coincide with the zeros of the canonical partition func- 
tion Z{B). For H = 0.1 mGauss only one distribution of ze- 
ros lying at negative temperatures is present (Fig. 0(a)) corre- 
sponding to the inverse occupation of the two isomers at neg- 
ative temperatures. With increasing magnetic field the shape 
of the distribution changes (Fig. ||(b)) and the influence of 
the poles of | (/i) | on the real temperature axis becomes vis- 
ible. While this effect is hardly seen for H = 0.1 mGauss, 
one is not able to distinguish the isomers due to their real- 
temperature values of (/i). For H — 1.2 mGauss and /? > 
the average magnetic moment equals 1/is, whereas we have 
(/i) — IQ^iB for negative temperatures. At higher mag- 
netic fields a second distribution of zeros corresponding to the 



structural transition between both isomers is seen. 

An inspection of both distributions of zeros reveals that the 
structural transition is of first order, i.e., for the classification 
parameters we have a = 7 = 0. Whereas for the transition lo- 
cated at negative temperatures we find a > and 7 7^ 0, i.e., a 
second-order transition, unless the magnetic field completely 
disappears. 

Another interesting example for first-order transitions, 
which become second order with growing size have been 
found by Proykova et al. [ETfl. There, two transitions have 
been found in TeFg clusters. They conclude that one of them 
is supposed to be of first-order, which becomes continuous for 
larger clusters, because they find two minima in the free en- 
ergy as a function of an order parameter which merge to one 
mimimum for larger cluster sizes. The arbitrary choice of the 
order parameter, however, might lead to different results p^]. 

Figure g displays the the magnetic susceptibility x within 
the complex temperature plane. The susceptibility is plotted 
versus the complex inverse temperature. The blurred peaks 
of X in Fig- |l(d) can be clearly related to "radiations" of the 
zeros onto the real axis. 

An inspection of the whole complex temperature plane re- 
veals the origin of the two maxima of the susceptibility at 
H = 0.5 mGauss on the real axis (cmp. Fig. ^(d)). For 
H — 0.5 mGauss the hump in x located at T « 3 K 
(/3 ~ 3800 1/eV) has its origin in the distribution of zeros 
lying at negative temperatures. Also the hump at T = 80 K, 
where earlier calculations have suggested that it is also related 



to a structural transition [171, has its origin in this distribution 



of zeros. For 2.0 mGauss the distribution of zeros lying at pos- 
itive temperatures contributes most to x for real temperatures 
at /? « 400 1/eV = T w 30 K corresponding to the structural 
transition seen in Fig. 0(d). 

Since the distribution of poles of the considered thermody- 
namic function (the distribution of zeros of Z{B)) is discrete 
the influences on the real axis might be shielded by zeros of 
these function. For example, the distributions of poles of (/i) 
and X ^"6 surrounded by distributions of zeros which are dif- 
ferent for (/i) and x- These thermodynamic functions are an- 
alytic everywhere except at the poles and zeros. Thus, if there 
is a zero in the vicinity of a pole near the real axis the decrease 
of (ii) or X cannot be compensated. 

In conclusion we have shown that the use of the complex 
inverse temperature plane has advantages to common investi- 
gations of thermodynamic functions. Within our classification 
scheme the order of a transition can be clearly identified. By 
means of a simple two-isomer model for small magnetic clus- 



ters we are able to identify two different types of phase tran- 
sitions. Furthermore, we found that signals in the magnetic 
susceptibility at positive temperatures might have their origin 
at negative temperatures. This also indicates that the inverse 
temperature is analytic at (3 — 1/T — and therefore should 
be used in calculations of thermodynamic properties. 

We thank E. R. Hilf for fruitful discussions and valuable 
comments. 
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